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Abstract 

The constraint distribution in non-holonomic mechanics has a dou- 
ble role. On one hand, it is a kinematic constraint, that is, it is a 
restriction on the motion itself. On the other hand, it is also a re- 
striction on the allowed variations when using D'Alembert's Principle 
to derive the equations of motion. We will show that many systems 
of physical interest where D'Alembert's Principle does not apply can 
be conveniently modeled within the general idea of the Principle of 
Virtual Work by the introduction of both kinematic constraints and 
variational constraints as being independent entities. This includes, 
for example, elastic rolling bodies and pneumatic tires. Also, D'A- 
lembert's Principle and Chetaev's Principle fall into this scheme. We 
emphasize the geometric point of view, avoiding the use of local coor- 
dinates, which is the appropriate setting for dealing with questions of 
global nature, like reduction. 
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1 Introduction 



Non-holonomic Mechanics. The universal formalism created by La- 
grange is not appropriate to derive the equations of motion for systems 
with rolling constraints, that is, this motion is not described by classical 
variational principles. Several systems with rolling constraints, like the ide- 
alized rigid ball rolling on a plane with only one point of contact and many 
others, are successfully described geometrically by distributions on config- 
uration space and the corresponding equations of motion are derived by 
D'Alembert's Principle, which has been the purpose of extensive research 
[s3 El HH El El EU EHi for more than a century (see also for instance 
|24l 121 171 IH] for a list of references and historical remarks). However, as we 
will see in the examples studied in the present work the dynamics of elas- 
tic rolling bodies is not generally described by D'Alembert's principle, even 
in those cases where the restriction on the motion is given by linear con- 
straints. On the other hand, second order constraints, that is, subsets 
of the second order tangent bundle rather than the tangent bundle of the 
configuration space, naturally appear in several examples. The purpose of 
the present work is to establish the basic geometric definitions and proce- 
dures within the general idea of the Principle of Virtual Work, generalizing 
D'Alembert's Principle to deal with nonlinear and higher order constraints. 
One of our main examples will be elastic rolling bodies, like pneumatic tires, 
where some second order constraints appear naturally. 

In D'Alembert's principle the constraint distribution has a double role. 
On one hand, it is a kinematic constraint, that is, it is a restriction on 
the motion itself. On the other hand, it is, in addition, a variational 
constraint. This perspective was already adopted in JH! where a general 
approach to non-holonomic constrained systems considered as implicit differ- 
ential equations was considered. There it was discussed that the kinematical 
constraints defining a submanifold on the tangent space of the configuration 
space of the system and the reaction or control forces described by a sub- 
bundle of the cotangent bundle of the configuration space, were independent 
entities and a condition for the compatibility of both ingredients was ob- 
tained. In this paper we will push forward this point of view by considering 
nonlinear higher order non-holonomic constraints, not only constraints on 
the velocities but on higher order derivatives. 

We will show that many systems of physical interest where D'Alembert's 
Principle does not apply, can be conveniently modeled by a Principle based 
in the introduction of both higher order kinematic constraints and higher 
order variational constraints as being independent entities. This in- 
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eludes, for example, elastic rolling bodies and pneumatic tires. Also, D'Alem- 
bert's Principle and Chetaev's Principle fall into this scheme. 

Our point of view is geometric, avoiding the use of local coordinates, 
which is appropriate for dealing global problems, like reduction. We also 
write equations of motion for systems with higher order constraints in an 
intrinsic fashion, using the natural structures of the tangent bundle and 
higher order bundles. 

Basic Notation As usual we will consider that the configuration space 
of a Lagrangian system is a smooth manifold Q of dimension n with local 
coordinates q l . We shall introduce higher order tangent bundles in order 
to deal with higher order constraints. Thus, by definition, two given curves 
in Q, say, qi(t) and q2{t), t G (—a, a), have a contact of order k at go = 
gi(0) = 52(0) if there is a local chart (99, U) such that gj(0) G U, for i = 1, 2, 
and Dl {tp o q{) (0) = Df [ip o q?) (0), for s = 0, k. This is a well defined 
equivalence relation, and the equivalence class of a given curve q(t) is denoted 
[q]( k \ For each q G Q, let T^ ] Q be the set of all [q]W such that q(0) = q , 
and let T^Q be the collection of all Tgg Q, for go G Q. It is well known 
(see for instance |19| . [§| and references therein) that T k : T^Q — > Q, where 
T k ([g]( fc )) = g(0), is a fiber bundle, called the tangent bundle of order k of 
Q. There are natural maps t' ' ' : 

T( fc )Q -> r«Q, for jfe,| = 1,2,..., given 
by r^' fc ) Qg]W) = [qf\ It is easy to see that T&Q = TQ. Also, we can 
identify T<®Q = Q, via [q](°) = g(0). 

In local coordinates, we have q = (q 1 , g n ), and, for s = 1,2,.., we 
denote q^ = (a 1 * 3 ', ...,q n '^) , where 

,<■<" = ^(0), 

where i = 1, n. Then we have, [q]^ = (g^, q^) ■ 

Denote by j k : T^Q -> T^'^Q) the canonical immersion defined by 
j k ([q]W) = [g( fc - 1 )]( 1 ) where g^" 1 ) is the lift to T^'^Q of the curve q, that 
is, the curve g^ 1 ) : (-a, a) -» is defined as = [gt] (fe_1) 

where gt(s) = g(i + s). 

In this paper, it will be useful to introduce, geometrically, the concept of 
implicit differential equations. This concept has often received less attention 
than the notion of an explicit differential equation in the differential geome- 
try literature (see |2S1 ^] ) • Geometrically, a system of implicit fcth-order 
differential equations is a submanifold M of T^Q and a curve 7 : I — ► Q is 
a solution to the differential equation M, if its fc-lift y( k '(s) G M for all s G /. 
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The implicit differential equation will be said to be integrable at a point if 
there exists a solution 7 such that its /c-lift passes through it. The integrable 
part of M is the subset of all integrable points of M. The system is said to 
be integrable if its integrable part coincides with M . A notorious algorithm 
has been developed to extract the integrable part of an arbitrary implicit 
differential equation [23], but it will not be the objective of this paper to 
discuss this issue for systems with higher order non-holonomic constraints 
and we will restrict ourselves to the description of the corresponding implicit 
differential equation, leaving the questions of the existence and uniqueness 
of its solutions for future discussion. 

In section |21 we describe a first example of the elastic rolling ball, where 
some of the features of the general procedure already appear. In sections |3] 
and^lwe show how to study Rocard's theory and also Greidanus's theory of 
a pneumatic tire (see ^3 EH1 H3 and also [21] ) as a non-holonomic system 
with higher order constraints and, motivated by the previous examples, in 
section [I] we establish a general principle for dealing with systems involving 
higher order constraints. The distinction between kinematic constraints and 
variational constraints as being independent entities is a key point to this 
discussion. In Section[5]intrinsic equations of motion for systems with higher 
order constraints are derived. In Section 6 further examples are provided 
and some basic results about reduction and the equations for Lagrangian 
systems with symmetries with higher-order non-holonomic constraints are 
discussed. 

2 A Simple Example: the Elastic Homogeneous 
Rolling Ball 

The main purpose of this section is to show an example that can be treated 
using D'Alembert's Principle and also using some other procedures involving 
different types of constraints, including second order nonlinear constraints. 
All those procedures are equivalent in the sense that they give equivalent 
systems of equations. 

Let us consider an elastic ball subjected to gravity and rolling on a plane. 
Without loss of generality we will assume that the radius of the ball is 1, 
for simplicity. For a static ball the contact between the ball and the plane 
is a circle, whose diameter was calculated by Hertz ^2j, see also ^H], page 
27. The effect of internal viscosity, adhesion and other dissipative forces 
is important in some cases [I], however, in the present section we shall as- 
sume that heat dissipation is small, in other words, we will consider only 
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the idealized model of a perfect elastic ball. Also, we shall consider only 
the important case where the circle of contact is small and the motion is 
quasistatic, which, in particular, implies that the zone of contact is approx- 
imately a circle of the same size as the contact circle in the static case (see 
jltjj). This also implies that the size and inertia of the flattened part of the 
sphere is negligible. Now we shall define the non-sliding condition. It is 
given by the condition that the points of the sphere belonging to the circle of 
contact cannot slide against the plane. It is clear that this has to be under- 
stood in an approximate sense since the exact solution of elasticity equations 
is not known in general, not even under the quasistatic assumption. More 
precisely, we accept the following approximate model. We assume that for 
all kinematic and dynamical purposes the ball is rigid, it has only one point 
of contact a with the plane, representing the center of the circle of contact, 
which does not slides, and the spatial angular velocity and the translation 
velocity combine in such a way that the velocity of points z of the surface 
of the ball near a have a velocity which is of order \a — z\ 2 . This is a rig- 
orous way of defining the constraint given by the non-sliding condition, in 
the case where there is a circle rather than a point of contact. It is easy 
to prove that, in fact, the non-sliding condition is satisfied if and only if 
the vertical component of the spatial angular velocity is 0, that is, U3 = 0. 
We emphasize that this model is realistic only for slow motion and small 
deformation. In agreement with all these physical assumptions we have the 
following geometric model. 

Kinematics of the Elastic Rolling Ball. The manifold Q = SO {3) x R 2 
is the configuration space for the model. A position of the system is given 
by a point (A, a) G Q, where a is the point of contact of the sphere with 
the plane representing in the approximation described above the center of 
the circle of contact. Let V = a be the translation velocity of the ball 
and let to = AA -1 be the spatial angular velocity uj = (u;i, u>2, W3), after the 
identification of so(3) with R 3 . We have V = a and u = AA~ l -AA^AA" 1 . 
The following two equations describe the non-sliding constraint 



The first equation represents the usual non-sliding condition for a rigid 
rolling ball while the second expresses the fact that there is really a circle 
of contact rather than a point, and that the points of that circle belonging 
to the sphere have zero velocity with respect to the plane, at least to first 





uj 3 = 0. 
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order approximation. The previous equations define a distribution, which 
is the kinematic constraint for the system of the elastic rolling ball. We 
will show that, provided that we accept higher order constraints, there are 
other equivalent ways of choosing the constraints all of them giving equiva- 
lent equations of motion. For instance, let the curve a(t) in the plane have 
curvature radius r(t). Then we define the constraint 



whose physical meaning is that the instantaneous motion of the sphere is a 
superposition of a rotation about some vertical axis, with angular velocity 
u>3, and the motion of rolling on the plane with speed 



and the point of contact is located at a distance r from the vertical axis. This 
is an example of a second order constraint, it is a kinematic constraint in 
the terminology introduced in section |I] and it is equivalent to the constraint 
(|2|). in the sense that it gives equivalent equations of motion, as we will 
explain later. However, as we have said before the non-sliding condition is 
satisfied only if r = oo, which of course implies = 0, or if u = 0. Equation 
(^Q) has the following consequence 



Let t and n be the tangent and normal vectors to the curve a(t). We have 



22 2,2 
r uj 3 = u{ +uj 2 , 



(3) 




(4) 




|V|n = ±(loi,lo 2 ) 



and also 



V = 



d\V 



V\ 2 



dt 



t + 



n. 



r 



Then we can deduce 




(5) 



VI s 



(6) 



r 



from which we obtain the constraint (j3J) in the form 




(7) 
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where the choice of the sign ± is the only one consistent with the standard 
choice for the direction of the normal n and the sign of 0J3 for the given 
physical description. We have a subset C C T^Q, given by (0) and (J7J), 
rewritten in terms of a, A, A, and A. This is a second order kinematic 
constraint. Observe that, in this example, the projection Tq' 2 ^ : T^Q — > 

TQ defines a distribution D C TQ, by D = Tq ' (C), which is given by (pQ), 
and that rewritten in terms of A, A, a and d, gives an expression linear in 
A and d. 



Dynamics of the Elastic Rolling Ball. The Lagrangian is given by the 
kinetic energy 

L(A, a, A, a) = \l{AA- l f + \m{o)\ 

where / is the moment of inertia of the ball with respect to any of its 
symmetry axis, and M is the mass of the ball. The dynamics of the elastic 
rolling ball is given by the following variational description, as we will see 
later, 

5 jT (^-I{AA- l f + lM(d) 2 )dt = (8) 

(5A(ti), 5a(ti)) = 0, for i = 0, 1 (9) 
(5A(t),6a(t)) ED {A(t)Mt)) , for all t (10) 

(i(t),a(t))e% )ia(t)) , for all t (11) 

co 3 = 0. (12) 

We will show that we can replace the last equation by equation (0 and we 
will obtain an equivalent system. We note that in this formulation the con- 
straints on the variations are the same as in the case of the rigid rolling 
ball (see for instance |24[ |2j). However, the kinematic constraints are 
not, in other words, the motion is effectively constrained by our choice of the 
last equation, namely, either equation (j2j) or equation 0. For any of those 
choices, we derive from the previous Principle a differential- algebraic 
system of equations and we will have existence and uniqueness of solution 
for those initial conditions compatible with the constraints. 

By applying the usual integration by parts argument, we obtain the 
equations of motion. However, as it already happens in the case of the rigid 
body, this is not completely trivial unless one is willing to use reduction 
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arguments, (see for instance jS] and [7]). We will postpone the details of the 
computation until Section 6. We obtain, 



(I + M)ui = 
(I + M)Co 2 = 
(I + M)lo 3 = 
(u 2 , -ux) = V 



(13) 
(14) 
(15) 
(16) 
(17) 



w 3 = 0. 



Of course this system is over determined, but it is correct. The fifth equation, 
which coincides with equation PJ). may be replaced by equation ((7J) and we 
obtain a system which is clearly equivalent. The first four equations are 
exactly the equations for the rigid rolling ball and they imply that u = 
and also that the translation velocity V is constant. We can show that there 
is solution provided that the initial condition {ujq, Vq) satisfies the constraints 
given by the last two equations and that this solution is unique. 

We must remark at this point that the only guiding idea to establish the 
previous procedure is the Principle of Virtual Work, and one should check 
that the final equations are consistent with the basic laws of mechanics, es- 
sentially Newton's Law, so the force should be equal to the rate of change 
of linear momentum and the torque should be equal to the rate of change of 
angular momentum. In the case of the elastic rolling ball the forces of the 
constraint must satisfy the following conditions: the resultant force exerted 
by the plane on the ball has a positive component in the vertical upwards 
direction while the torque has a zero horizontal component. All this is obvi- 
ously compatible with the previous system of equations. Moreover, the same 
equations can be derived by an elementary exercise in rational mechanics. 
We observe that preservation of energy is satisfied in this example. As a 
final remark to this example we observe that even if the constraints (fljl. (J2J 
are linear, we have not applied D'Alembert's Principle. However, it will be- 
come clear at the end in section |H1 that D'Alembert's Principle gives correct 
equations of motion in this example, and it is perhaps the best procedure 
in this case since it produces a non-overdetermined system. Showing that 
it is not always the case that D'Alembert's Principle can be applied is part 
of the purpose of the present work. It is also clear from what we have ex- 
plained so far that, for a given system, there is in principle the possibility of 
introducing several classes of higher order constraints which are equivalent 
in the sense that they lead to equivalent equations of motion. 

The case of the nonhomogeneous elastic ball and also the case of the 
nonhomogeneous viscoelastic ball could be interesting, for instance because 
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of possible applications to spherical robots, and can be treated with the 
methods of the present work. In particular, the non-sliding condition @ 
will be part of the kinematic constraints. The case of the symmetric elastic 
or viscoelastic rolling ball, in which two of the three moments of inertia of 
the ball are equal, presents an extra symmetry and we can expect that some 
kind of reduction by this symmetry will help to understand the behavior of 
the reduced variables such as the angular momentum. The case of the rigid 
symmetric rolling ball has been studied in [5]. 

3 An Example of Nonlinear Higher Order Non- 
holonomic Constraints 

In the example of the elastic rolling ball described in the previous section the 
second order constraint gives rise to a distribution D defined by (^Q) which 
provides a restriction for the variations to obtain some of the equations of 
motion. The rest of the equations of motion are the ones given by the same 
distribution, plus an extra equation provided by the nonlinear second order 
constraint or, equivalently, by the linear constraint (J2J). This gives a pro- 
cedure whose correctness in the example under consideration is established 
by the fundamental principles of mechanics. 

Rocard's Theory of a Pneumatic Tire. Before we try to establish 
any general procedure we will describe another example where the restric- 
tions, both kinematic restrictions and restrictions on the variations, are of 
an entirely different nature. This is the simplified model of a pneumatic tire 
rolling on a plane according to Rocard's theory, as described for instance 
in |27j . |26j . |24j . For simplicity we shall study the case of a single elastic 
pneumatic tire whose plane is constrained to remain vertical while it rolls 
without sliding. The zone of contact of the pneumatic tire with the plane 
is a small surface with a central point of contact x = (xi,X2), which for 
simplicity we will assume that it coincides with the projection of the center 
of the wheel on the plane. The non-sliding condition means that the velocity 
of the points of the tire belonging to the zone of contact with respect to the 
plane is zero. In an approximate sense this non-sliding condition implies 
that the vertical component of the angular velocity of the small piece of sur- 
face of the pneumatic in contact with the floor is zero. However, contrary to 
what we have assumed for the homogeneous elastic rolling ball, the fact that 
the vertical component of the angular velocity of the zone of contact is zero 
does not mean that the vertical component of the angular velocity of the 
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plane of the tire is zero. This is because according to Rocard's theory the 
elasticity of the material allows for a small angle e between the axis of the 
zone of contact(an oblong-like symmetric zone), which is assumed to have 
the direction of x, and the plane of the wheel. We will call K the corre- 
sponding constant of elasticity. It turns out that the non-sliding condition 
for the small zone of contact is not the relevant constraint. Instead, there 
will appear another second order constraint of a different nature. Finally, 
we must remark that the previous description of Rocard's theory gives only 
an approximation, and for more accurate results one must have into account 
some other observed effects. For instance, the projection x of the center of 
the wheel onto the plane is not exactly the center of the zone of contact, 
which produce a small torque not taken into account in the simplified model 
described above. Part, but only part, of this problem is taken into account 
in the simplified version of Greidanu's theory described later in the present 
work. 

Taking into account all the physical considerations explained above we 
will describe Rocard's theory by the following geometric model. For all 
kinematic and dynamical purposes the wheel is simply an undeformable 
disk kept vertical and rolling on a plane, where the point of contact is x = 
(x±,X2)- We choose once for all a normal vector N = (— sin 6*, cos^) rigidly 
fixed to the wheel. Then the angle between the plane of the wheel and the 
x\ axis is 9. The angle between the velocity vector x and the plane of the 
wheel is called e, with the physical meaning that we have explained before. 
Therefore, the angle between the axis x\ and x is — e, and the vector 
n, normal to the trajectory of the point x and pointing in the direction 
of the concavity of the curve, is n = (— sin(# — e), cos(# — e)) . The angle of 
rotation of the wheel about its own axis is called ijj. In order to obtain precise 
formulas one should be careful about the sign conventions. Positive angles 
in the x\X2 plane satisfy the usual convention. Thus the angle between the 
x\ axis and the X2 axis is, by definition, (l/2)w while the angle between the 
X2 axis and the x\ axis is — (l/2)-7r. The sign for the angle ip is established 
by the convention that the vector angular velocity is of the form ipN. The 
configuration space of the system is Q = T 3 x B 2 , and a generic point is 
Q = {Qi,Q2, <?3, <?4, 9s) = (ip, Q, e, xi,x 2 ). The Lagrangian is given by 

L(q, q) = \H 2 + \ J9 2 + l -Mx 2 - ^Ke 2 , 

where / is the moment of inertia of the wheel with respect to its axis, J is the 
moment of inertia of the wheel with respect to any one of its diameters, M is 
the mass of the wheel and K is the constant of elasticity introduced before, 
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which by definition satisfies T = —Ke, where T is the vertical torque. The 
kinetic energy due to the velocity of rotation e of the small flattened piece 
of material about the zone of contact is small and we will assume that it is 
for simplicity, which is also in agreement with general standard assumptions 
for this kind of approximate models, |24| . 

Next we shall describe the kinematic constraints and the variational 
constraints. The kinematic constraint Ck, is given by the equations 

xx = ^>cos(0-e) (18) 
x 2 = ip sin(0 - e) (19) 
-^tge + ^{6-e) = (signVO-^rtge. (20) 

The first two equations represent the non-sliding condition for the center of 
the zone of contact, and they are the same as the ones that appear in the 
case of a rigid rolling disk, or wheel, except for the small angle e. We should 
emphasize that here we are working to first order approximation only, which 
means that powers of e greater than 1 may be neglected. The last equation 
comes from Rocard's condition, 

\F\ = asm. |e|, 

where a is a positive physical constant and F is the force normal to the 
wheel exerted by the floor, while the wheel is rolling with nonzero velocity. 
More precisely, F is the N component of the centripetal force, that is we 
have F =< Mx, N > . The sign conventions are encoded in the following 
more precise version of Rocard's formula 

F = (sign ip)a sine, 

where e must be interpreted as being the angle between the normal n to 
the curve and N if F > while it must be interpreted as being the angle 
between n and —N if F < 0. Recall that Rocards's formula is valid for 
e close to only. A couple of remarks is in order for future use. First, 
as we have said before, Rocard's theory is valid modulo infinitesimals of 
order (sine) 2 . Second, with the previous sign conventions and according to 
Rocard's formula it is not difficult to show that e{9 — e) > 0. It also follows 
from the expression of Rocard's formula given by (|2U|) that for e = the 
curve x(t) must have a point of inflection, that is 6 — e = 0. 

It is clear that (|2()|) involves the first and second derivatives of some of 
the variables with respect to time, moreover, the dependence on the first 
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derivatives is nonlinear, therefore it is far from the typical constraints of 
D'Alembert type. To obtain equation Q2U|) we may assume, without loss of 
generality, that ip > 0. We simply differentiate (JT8|) and (|T§)) with respect 
to time, and replace in the equation (sign ip)a sine =< Mx,N > . Now 
let us consider the following variational constraints Cy, to be imposed on 
variations 5q 

Sip cos 0-Sx 1 = O (21) 
5ip sin 9-5x2 = (22) 
59 - 5e = 0. (23) 

Consider the curves q(t) satisfying 

5 [ L(q,q)dt = 0, 

Jto 

for variations 5q satisfying 5q(ti) = 0, for i = 1,2, and also the variational 
constraints Cy . Those curves are the ones satisfying the following dynamic 
equations 

I-ij) + Mil cos 9 + Mx 2 sin 9 = (24) 
JO + Ke = 0, (25) 

obtained by the usual integration by parts arguments. These dynamic equa- 
tions give balance between forces of the constraint and rate of change of 
momentum. The resultant of the forces exerted by the plane of contact on 
the wheel has positive upwards vertical component which is compensated 
by gravity, while the horizontal component, which is given by Mi, is de- 
composed in the directions (cos 9, sin 9) and (— sin#, cos 9). The first one 
is compensated by the rate of change of the angular momentum lip and 
the second is compensated by the non-sliding constraint force. The vertical 
component of the torque of the forces exerted by the plane on the wheel 
is Ke which is compensated by J9. The other components of the torque 
are automatically compensated because we are assuming that the wheel is 
forced to remain vertical. The system of dynamic equations together with 
the kinematic constraints equations Ck completely describe the motion of 
the wheel. 

In the previous example, we should emphasize, again, the distinction 
between kinematic constraints and variational constraints. They are 
conceptually different, and this difference is implicit in the usual statement 
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of the Principle of Virtual Work. However, in the literature this distinc- 
tion is usually not emphasized, and for good reason, since in those cases 
where D'Alembert's principle can be applied the variational constraints and 
the kinematic constraints coincide. Non-holonomic systems that cannot 
be treated using D'Alembert method have been considered for instance by 
Chetaev [IJJj where a procedure to deal with general first order nonlinear 
constraints is devised (see also EH])- In Marie [22] it is clearly stated 
that constraint forces cannot be derived in general from the kinematic con- 
straints and have to be added as part of the physical description of the 
system. Furthermore in |13| it was explicitly stated a formulation for first 
order Lagrangian and Poisson nonholonomic systems where kinematic con- 
straints and constraint forces are given as independent entities. 

In the case of the elastic rolling ball the forces of the constraint are 
normal to the direction of the motion of the ball and there is no dissipation 
of energy. However, for a viscoelastic rolling ball there is certainly dissipation 
of energy and the component of the force of the constraint in the direction of 
the motion can be calculated using results from 0J. This kind of system can 
also be approached using the kind of generalization of D'Alembert's principle 
described in section 0J The rate of dissipation of energy for a pneumatic 
tire rolling according to Rocard's theory can be easily calculated. Since the 
energy is given by E = (l/2)Iip 2 +(l/2)J9 2 +(l/2)Mx 2 +(l/2)Ke 2 , using the 
kinematic constraints (|18|). (|19j) and the dynamic equations derived before 

we can show after some easy calculations that E = — (^Mip 2 + K^j e{9 — e), 

modulo infinitesimals of order e 2 . Since e(9 — e) > as we have explained 
before we have E ^ 0, which means that in general there is dissipation of 
energy. The limit case e = gives E = 0, which reveals that Rocard's theory 
does not take into account the relatively small dissipation of energy that 
occurs when the tire rolls in a straight line. To prove the previous formula we 
proceed as follows. We can easily see that E = Iip'ip + J99 + Mx-x + Ke€. By 
differentiating ()18|) and ()19|) we can easily see that x ■ x = ipip and from this 
and the dynamic equation (|25j) we obtain (I + M)ipip — Ke(9 — e) = 0. Using 
(|18(1 . ((T§|) and we obtain, modulo higher order infinitesimals, that (/ + 
M)ip = -Mipe(9 - e) therefore (I + M^tp = -Mip 2 e(9 - e), from which we 
finally obtain E = — ^Mip 2 + K^j e(9 — e). ^From a general point of view we 

may say that the distinction between variational and kinematic constraints 
implies that the infinitesimal work of the constraint forces in general does not 
vanish for some admissible infinitesimal displacements, which is the reason 
why the forces of the constraint may produce work. 

In the next section we define a class of non-holonomic systems with 
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higher order nonlinear constraints based on the introduction of both kine- 
matic and variational constraints. We will also show that procedures like 
D'Alembert's Principle or Chetaev's procedure fall into this scheme. We pro- 
pose that questions of a general nature on non-holonomic systems, like reduc- 
tion by the symmetry, Legendre transformation, and many others should be 
approached for the general case of higher order constraints using the scheme 
based on the introduction of both kinematic and variational constraints. 

4 A Principle of Virtual Work for Lagrangian Sys- 
tems with Nonlinear Higher order Non-holonomic 
Constraints 

Let Q be a configuration space of dimension n and let L : TQ — > R be a given 
Lagrangian. Then we have the Euler-Lagrange operator ££ : T&Q -> T*Q 
which is given in coordinates by 

^*l M ^=(li(W M )-|(H w ))««'- 

A kinematic constraint of order k is, by definition, a subset Ck Q 
T<MQ, for some k = 0, 1, 2, ... The subset Ck is often defined by equations 
R K ({ q ]( k )) = 0, where R K : T&iQ -» W, for some r = 1,2, .... For example 
if k = and Rk is a submersion then Ck is a nonsingular holonomic con- 
straint. If k = 1 and RK(q,q) = RKi(q)q l defines a distribution of constant 
rank, we have the typical situation of D'Alembert's Principle. If Rk(q, q) is 
a general function we have the situation considered by Chetaev ^U] . In the 
case of the elastic rolling ball we have, if we choose the constraint given by 
equation (J2J) as we have explained before, n = 5, k = 1, r = 3, and 

Rk(A, a, A, a) = (o>2 — &i, — a 2 , 0J3). 

Alternatively, as we have explained before, if we choose the constraint given 
by equation (J7J), we have, n = 5, k = 2, r = 3, 

Rk(A, a, A, a, A, d) = (U2 — di, — oji - 02, ^1^2 — — u)$(oj\ + 

In the case of the Rocard's theory of a pneumatic tire, we have n = 5, k = 2, 
r = 3, and 

R K (^p,0,e,xi,x 2 ,ip,0,e,xi,X2^ (26) 

= [x\ — ificos(9 — e), X2 — Tpsm(9 — e), — V'tge + tp(9 — e) — (sign^)— tge 

(27) 
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A constraint on the variations of order I is a subset Cy C T®Q xq 
TQ denned by equations Ry {\q\^\ Sqj = where Ry is linear in the vari- 
able 5q, so we shall write as usual Ry ([q]^ l \ Sq) = Ry ([q]^) -Sq or, in 
coordinates, Ry ( [<?] (/) , <5<?) = Ryi ([q] {l) ) ■ Sq\ For each G T^Q, we let 
C v {[q]W)={6q€TQ:([q]W,6q)€C v }. 

Statement of the Principle. The main object defined in this paper is 
the class of Lagrangian non-holonomic systems defined by data (L, Ck, Cy) 
whose dynamical equations are derived by using the variational principle 

5 f L(q,q)dt = 0, 

where variations 5q are restricted by 5q G Cy ([<j»]®), or, equivalently, 
Ry ([q]^) ■ Sq = 0. Then the equations of motion are given by the dy- 
namical equations 

e£i([q]V)) € R V ([q]®)° 
and the kinematic constraint equations [q]® G Ck or, equivalently, 

R K ([q] (k) )=0. 

Equations of motion will be derived in the next section. 

The previous Principle, which is contained in the general idea of the 
Principle of Virtual Work, imposses, through the dynamical equations, re- 
strictions on the forces of the constraints. But, contrary to what happens 
with D'Alembert's Principle, the forces of the constraints derived from the 
Principle stated above will in general produce work. 

The class of higher order non-holonomic systems just defined contains 
several important classes of non-holonomic systems. For example, for the 
class of non-holonomic systems that are tractable using D'Alembert's princi- 
ple we have, by definition, k = 1, I = and Ck is the distribution where for 
each q G Q the space of the distribution is Cy(q) C TQ. Thus, the kinematic 
constraint and the constraint on the variations essentially coincide in this 
case. In the case of nonlinear kinematic constraints considered by Chetaev 
given by Rx(q,q) = we have I = 1 and the variational constraints are 
defined, according to Chetaev, by 

d i -\ x dR K {q,q) x 

Rv(q,q) ■ oq = — ^ oq. 
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Remark 4.1 In the mathematical literature one finds some examples of 
higher order constraints in non-holonomic problems (for instance see |281 1251 
1151 I29| ) . In the previous references an extension of the Chetaev principle for 
kinematic second order constraints is applied, namely, 

(R K )i(q, q, q) = 0, 1 < i < m 

and variational constraints Ry are derived from the kinematic constraints 
by 

Rv{q, q, q) ■ $q = —q~ - -6q = o 

In the case of the elastic rolling ball the variational constraints are given 
by IjlOJ) . In the case of the pneumatic tire according to Rocard's theory 
the kinematic constraints are given by (jlHjl . (|19j) . (1201) and the variational 
constraints are given by (JUJ, (JUJ), (|25|) . 

We emphasize once again that the notions of kinematic constraints and 
variational constraints are independent and one should not attempt, for 
instance, to derive variational constraints from kinematic constraints by a 
universal procedure. In order to illustrate further the necessity of such 
a point of view we will describe next the example of Greidanus's theory 
of a pneumatic tire, where the kinematic constraint defines a distribution 
like in D'Alembert's Principle but the variational constraints are not given 
by the same distribution, therefore they are not the ones prescribed by 
D'Alembert's Principle. 

Pneumatic tires according to Greidanus Several approaches to the 
dynamics of a pneumatic tire like those of Rocard, Greidanus, Keldys and 
others can be found in |27j . |26j . |24j . To describe Greidanu's approach 
we shall consider the simpler setting of Rocard's approach described before, 
but this time we allow, in addition, for a lateral deformation £. The absolute 
value of the quantity £ is the distance between the projection of the center 
of the wheel on the plane (x\,xz) and the center of the zone of contact. In 
the Rocard's approach described above the value of £ is 0. We must remark 
that we are considering in this paper only the case of Greidanus's theory in 
which the wheel is kept vertical. The physical reason for the appearance of 
the displacement £ is of course the lateral deformation due to the centrifugal 
force given the elasticity of the material. 
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The kinematic constraints are 



x\ = ip cos(# — e) 
±2 = ip sin(# — e) 
9 -e = ip(a£ + (3e). 



(28) 
(29) 
(30) 



The first two equations are the same as in Rocard's approach. The last one 
expresses the fact that the curvature of the trajectory of the center of the 
contact zone is, for a given speed of rotation of the wheel, proportional to 
a linear combination of the deformation parameters £ and e, where a > 
and P > 0. This replaces Rocard's constraint. We see that the kinematic 
constraints define a distribution. The variational constraints are 



These variational constraints are different from the kinematic constraints, 
therefore we are not using here D'Alembert's Principle. The projection of 
the center of the wheel on the plane is the point (2/1,2/2) given by 



It is more convenient to calculate the kinematic constraints and the varia- 
tional constraints in terms of y\ and 2/2 instead of x\ and x 2 - The kinematic 
constraints are 



50-6e = 0. 



Sxi = 5ip cos 8 
5x2 = Sip sin 



(31) 
(32) 
(33) 



2/1 = x\ + ism.9 
2/2 = x 2 - £ cos 9. 



(34) 
(35) 




(36) 
(37) 
(38) 




9 



— e 



The variational constraints are 



S yi = 5il) cos 9 + 5£ sin 9 + £(cos 9)59 
5y 2 = 5tp sin 9 — 5^ cos 9 + £ (sin 9)59 
59-5e = 0. 



(39) 
(40) 
(41) 



The Lagrangian is 



+ l -M{(y 1 f + (y 2 ) 2 ) 



1 9 1 9 

r e ~ 2^ • 
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Then, equations of motion are given by kinematic constraints (|36[). (|37|). 
1)38(1 and dynamic equations 

7$ + M yi cos 9 + My 2 sin 6> = (42) 
J0 + M£yi cos 6» + sin 6> + /3e = (43) 

-Mj/i sin 6> + My 2 cos - a£ = 0. (44) 

We can easily check that the previous equations represent the balance be- 
tween rate of change of momentum and forces of the constraints. 

For high values of a the deformation £ remains small. Moreover, for a — > 
oo we have £ — > and the dynamic equations (|42l) . (|43|) of Greidanu's theory 
become the equations (f2~I|) . (f23|) of Rocard's theory, provided that K = (3. 
Using this and the fact that the two first kinematic constraints (|18|) . (|19|) of 
Rocard's theory coincide with the first two kinematic constraints (|28j) . (|29j) 
of Greidanu's theory and also the fact that for a — > oo the mechanical energy 
E for both theories tend to the same value, one can prove, proceeding as in 
the case of Rocard's theory, that at least for high values of a a pneumatic 
tire moving according to Greidanus theory is a dissipative system. This 
shows that D'Alembert's Principle does not provides a good model for this 
kind of system., even though the kinematic constraints are linear. 



5 Equations of motion 

Let us recall some basic facts of the geometry of the tangent bundle. The 
vertical endomorphism S is defined in local natural coordinates (q A , q A ) 
on TQ by 



(9) dr A 

The Liouville vector field A on TQ is locally defined by 

d 



A = q / 



dq' 



A second order differential equation is a vector field T on TQ such that 
S'(r) = A. We have the following local expression for T: 



■ a d a. .. d 



An integral curve of T is always the tangent prolongation of its projection 
q(t) on Q, called a solution of T. It satisfies the following explicit system 
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of second order differential equations: 



We also note that the kernel and image of S consist of vertical vector fields. 
Moreover, S acts by duality on forms and the kernel and image of S* consists 
of horizontal 1-forms. 

Given a lagrangian function L : TQ — > R, we construct the two-form 
col = —d(S*(dL)) on TQ, and the energy function El = AL — L (see HO]). 
A remarkable property of S and ujl is the following is^L = 0, or, in other 
words, 

S* o Cjl = — o S 1 , (45) 
where Cjl denotes the map T(TQ) — > T*(TQ) defined by contraction with 

LO L - 

Observe that if L is regular, then is a symplectic form, and there is 
a unique vector field Tl satisfying 

ir L u L = dE L , 

or, in other words, Tl is the Hamiltonian vector field with Hamiltonian 
energy El- It is well known that Tl is a second order differential equation 
on TQ, namely, the Euler-Lagrange equations for L. 

Without the regularity condition, the Euler-Lagrange equations form a 
system of second order differential equations in Q, in implicit form, that is, 
a submanifold D2 of T^Q, determined by: 

D 2 = {w G T^Q I i h(w) u L (r^ 2 Hw)) = dE L (r^(w))} (46) 

or, in other words, 

D 2 = {w G T (2) Q I £C{w) = 0} . 

The class of higher order non-holonomic systems studied in this paper, are 
determined by data (L,C^,Cy). Next we will show that the equations of 
motion of this kind of systems is a system of implicit /cth-order differential 
equations. In what follows, and without loss of generality, we will always 
suppose that k > I and k > 2. 

In our case the constraint on the variations are determined by a subset 
C v C T®Q xqTQ. Therefore for each point [q]® we obtain the annihilator 
C%r([q]W) C T*Q of C v ([q]®). Denote by F v ([q]®) the subspace of T*(TQ) 
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determined by F v ([q]®) = (r Q )*(C£([g]W)). Now, we shall define the subset 
of T^Q: 

My = {[ 9 ]« g T k Q | i M[q] i^ L ([q] {1) ) ~ dE L ([q] W ) £ F V ([?F)} ■ 

Therefore, the non-holonomic system associated to (L, C# , Cy), determines 
a fcth-order implicit system given by the submanifold Mkv = Ck fl My. 
The solutions of the problem (L, C/<-, Cy) are the curves 7 : / — > Q such 
that 7W C M K y. 

6 Further Results and Examples 

The scheme generalizing D'Alembert Principle, for the case of higher order 
constraints described in section|I]is not of course the most general case. It is 
not the purpose of the present work to expose the most general possible for- 
malism, but on the contrary, to provide a scheme which is useful in a variety 
of problems in mechanics. This scheme is also useful to deal with impor- 
tant questions of a general character in mechanics, like reduction, Legendre 
transformation and others. Some of these questions will be the purpose of 
future work and in this section we will consider some partial results only. 

Reduction of Invariant Systems with Higher Order Constraints on 
a Group. In this paragraph we explain how to reduce invariant Lagrangian 
systems with higher order non-holonomic constraints on a group. The more 
general case of systems on a principal bundle will be the purpose of a future 
work. However, in the present section we will show how to proceed in an 
example where the bundle is trivial, which illustrates some of the features 
of the general theory. Assume that the configuration space is a group G 
and that the Lagrangian L, the kinematic constraint Ck and the constraint 
on the variations Cy are left invariant. For right invariant systems we can 
proceed in a similar way. For each r = 1, 2, ... we have an identification 

a r : T^G/G -> rg, 

where rg = g ... q, is the direct sum of r copies of g. This identification 
is uniquely defined by the map [g]^ — > [v]^ r ', where v = g~ l g, and [v]^ = 
(i;( ), -i/ 1 ), ..."iA" 1 )) , where, by definition, we have, 

v {l> = — -v , 
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for r = 0, 1, ...r — 1. Under the identification a&, the quotient of the kinematic 
constraint Ck/G, becomes a subset, called reduced kinematic constraint, 
£k Q kg. Similarly, for each r = 1,2,... we have an identification 

P r : (t^Gx g TG) /G^rgeg, 

This identification is uniquely defined by the map {[g]^ r \Sg^ — > (H^^r/) 
with [f]( r ^ = (dC'^W, ...iA -1 )) , as before, and 77 = g~ 1 5g. Under the 
identification (3i, the quotient of the constraint on the variations Cy/G, be- 
comes a subset, called reduced variational constraints, €d Q Iq © 0- 
Since the equations Rx(\g\^) = and Ry([g]( l \5g) = are invariant, 
we have reduced equations !SHkt([v]^) = and 9fy ( [v ] ® , rj) = 0. Since 
RviM^^g) = R v • Sg is linear in Sg, we have that 9fy ([g]^) ■ n 

is also linear in 77. The Lagrangian L gives rise to a reduced Lagrangian 
I : q — ► M. We have the following theorem 

Theorem 6.1 The following conditions are equivalent 

(i) The curve g(t) satisfies 

5 I' L(g,g)dt = 0, 
J t 

for all 5g such that Sg(t) G Cy , /or a// i G [io>*i] f equiv- 

alently Ry ([g]^(t),Sg(t)) = /or a// f € [t , Sg(U) = 

/or i = 0,1; [g] (k) (t) G (equivalently R K ([g] {k) (t)) = /or a// 

t g [t ,*i];. 

(ii) The curve g(t) satisfies the equation 

/or a// <5</ stic/i i/taf <5^(i) G Cy > / or a ^ * e [*o>*i] (equiv- 

alently Ry ([g]V\t),5g(t)) = /or a// t G [i ,*i],) «^ Sgfc) = 
/or i = 0,1; [g]( k) (t) G Cx (equivalently R K ([g]^(t)) = /or all 

t€ [t ,h]). 

(Hi) The curve v(t) = 1 (i)<)(f ) satisfies 
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for all 5v = r] + [v,rj] where r)(t) G £y ([«]"'(<)) /or a// i G [*o>*i] 
(equivalently 5Hy (H^(t),7?(i)) = /or a// i G [*o,*i]^ ari ^ ?7(*i) = 0, 
for i = 0,1; [u] (fe) (t) G £ K (equivalently V\ K {[v] (k) (t)) = /or a// 

*€[t ,ti];. 

(ro^ T/ie curve v(t) = g (t)g(t) satisfies the equation 




for all 7] such that rj(t) G Cy (M^(*)) f or a ^ * ^ [*0i^l] (equivalently 
9V (M W ( i )> r /(*)) = / or a ^ G [*o,*i]J andriiU) = 0, /or i = 0,1; 
H (fc) (*) G £x (equivalently VK K (M (fe) (*)) = /or a// t G [t ,h].) 

The proof of this theorem can be performed proceeding as in [J]. The idea 
of the proof is simple. Given a curve g(t) such that [g]( k \t) G C'k for all 
t G [ioj^l] we take variations 5g(t) = g(t)r](t) for all t G [tcn^i] such that 
5g(t) G Cy ([g]®(t)) for all t G [*o,^l]- Since v(t) = gr 1 we can easily 

check that 5v(t) = r](t) + [v(t),rj(t)]. The rest of the proof follows by keeping 
track of the reduction of both the kinematic constraints and the variational 
constraints. 

Symmetry of the Elastic Rolling Ball. An interesting case occurs 
when, for each [g] W , Cy ( [g] ® ) depends only on g giving rise to a distribution 
D on G. This happens in the case of the rolling ball studied in section [21 Let 
us see how the previous theorem applies to this case. First of all we observe 
that the configuration space is the direct product group SO (3) x M?. Since 
we are assuming an homogeneous ball the kinetic energy Lagrangian is not 
only left invariant but also right invariant. This is important because the 
constraints are also right invariant. We can thus reduce by the right action 
of the group on itself. For r] = (a, w) and taking into account that the Lie 
bracket in so (3) is minus the standard one because we are reducing by right 
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actions, we have 




+ -MV 2 )dt = 
2 J 

[u),a] 

for i = 0, 1 



(47) 



a(U) = 0, 
5V = w 
w(U) = 0, 



5uj = a — 



w = (a 2 , 

U = (W2, 
^3 = 0. 



for % 




0,1 



(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
(54) 



Equations (|15|) , (j5U|) , (|5T|) represent the reduced variational constraints while 
equations (|52|). (|53j). (f51|) represent the reduced kinematic constraints (as 
we have explained before equation (|54j) can be replaced by u)2^i — oj\uj2 = 
iOs(u>f +lo 2 )). We obtain the equations of motions written in section[21 that is 
equations (|T3J), (JUJ), (CSJ), CU, OH)- The reduced version of D'Alembert's 
Principle consists of all the previous conditions plus the condition 03 = 
0, which of course corresponds to the kinematic constraint 0^3 = 0. The 
D'Alembert equations are (JT3J), (JHJ), CHI)' <[T7|). 

Rigid Ball Rolling on a Moving Plane. For dealing with examples 
where the configuration space is a principal bundle rather than a group 
and the constraints and also the Lagrangian are invariant we need to gen- 
eralize the previous theory, which we plan to do as part of future works. 
However, some simple examples can be worked out directly as we will see 
next. Let us consider a rigid ball rolling on a plane while this plane is 
being continuously deformed according to the law ipt : M 2 — > M 2 . The Eu- 
lerian velocity is vt(x) = <pt V^ix) an d we will assume that vt(x) = v(x) 
is independent of t. For a rigid ball rolling on a fixed plane, that is when 
v (x) = 0, the system is governed by the D'Alembert Principle which in 
this case is like the Principle of Virtual Work described in section for 
an elastic ball except that one should eliminate the kinematic constraint 
LJ3 = 0. When v (x) 7^ there is an extra force since the point a of the 
ball which is in contact with the plane, is moving with velocity v (a), that 
is, the kinematic constraint becomes (ui2,— <^i) = a — v(a). By differen- 
tiating with respect to t we obtain (cl^,— = a — Dv(a).a. Using this 
it can be easily seen that the force exerted by the floor on the ball is 
M ((^2, — ^1) + Dv(a) ■ (lj2, —oj%) + Dv(a) ■ v(a)) . Equations of motion can 
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be easily derived by direct application of the basic rules of mechanics and 
we obtain 



(J + M)(w2,-u;i) 



MDv(a) ■ [(w 2 ,-wi) + u(a)] 



(55) 
(56) 



UJ3 = 



Now we want to obtain the same equations using the formalism of the Prin- 
ciple stated in section |1J As in the case of the elastic rolling ball this is not 
straightforward, which emphasizes the advantages of having a way of reduc- 
ing by the symmetry as we will show next. The example under consideration 
is invariant with respect to the right action of SO(3) only because in this 
case the kinematic constraint is not necessarily invariant under translations. 
As we have said before in this simple example a general theory of reduction 
for systems on a principal bundle is not needed. Moreover, it is not difficult 
to prove directly that the following reduced Principle of Virtual Work gives 
the correct equations of motion 



Equations (|58|) . (|59|) and (|60j) represent the variational constraints while 
equation ()61|) is the kinematic constraint. 
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